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Abstract 

We introduce a notion of non-local almost minimal boundaries similar to that introduced by Alm- 
gren in geometric measure theory. Extending methods developed recently for non-local minimal 
boundaries we prove that flat non-local almost minimal boundaries are smooth. This can be viewed 
as a non-local version of the Almgren-De Giorgi-Tamanini regularity theory. The main result has 
several applications, among these C 1,a regularity for sets with prescribed non-local mean curvature 
in L p and regularity of solutions to non-local obstacle problems. 

Keywords: minimal surfaces; almost minimal surfaces; non-local elliptic equations; non-local 
mean curvature; geometric measure theory; integro-differential equations. 



1. Introduction 

In this article we introduce a notion of non-local almost minimal boundaries and study their reg- 
ularity near flat points as well as their singularities. In rough terms, a non-local minimal boundary 
corresponds to a set whose characteristic function minimizes a non-local energy functional, namely 
a fractional Sobolev norm (specifically H s / 2 , s < 1); they generalize classical minimal surfaces of 
codimension one, in particular, it has been observed that when s — 1 they are nothing but minimal 
surfaces. 

Non-local minimal surfaces appeared recently in the study of phase transitions with long range 
interactions, as equilibrium configurations. This was shown in 6], via a level set type scheme (see 
also Q). The regularity theory of non-local minimal surfaces was developed in Q, there it was 
shown that such surfaces are smooth except a singular set of at most dimension n — 2 (n is the 
space dimension), we follow closely the methods developed there. 

Our main result states that if an almost minimal boundary(for a non-local energy functional) is 
flat enough near a point, then in a neighborhood of that point it is smooth, i.e. it has a continuously 
changing normal. Namely: 

Theorem 1.1. Suppose E is ( J s , p, 5) -minimal in B\, where p satisfies assumptions Al through 
A3 (Section 2). There exists 5q = So(n,s,p) such that if E is So-flat in B\, then dE is C 1 in Bi/ 2 - 
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Moreover, we show that the set of points where dE is not C 1 is at most n — 2 dimensional. 

Theorem 1.2. If E is (J 8 , p,d) -minimal in O, then the Hausdorff dimension of the singular set 
Tie C dE HQ is at most n — 2. 

For the definition of a ( J s , p, 5o)-minimal set and <5-flatness, see Section 2. 

The classical theory of almost minimal boundaries says that a set E whose perimeter does not 
decrease too much under small perturbations has a smooth boundary, except perhaps for a lower 
dimensional set of singularities. This is known as the Almgren-De Giorgi-Tamanini regularity the- 
ory (see Chapter 4 of [H or [lj|, also the classical reference Q). The term "almost" refers to the 
fact that the energy does not decrease too much under small perturbations. In the present paper we 
prove that a similar statement holds if instead of the perimeter we consider the energy associated 
to H s / 2 , which is a non-local functional. It is worth mentioning also that the regularity theory for 
fully non-linear elliptic equations has been developed recently for non-local operators [3] . 

The proof of the above theorem follows substantially the techniques in Q which are based on 
Savin's proof of regularity of minimal surfaces and level sets in phase transitions (see 13 1 and [3]), 
they are also closely related to work of Caffarelli and Cordoba 

To prove regularity of non-local minimal surfaces there are three basic steps [H, Hj]: 1) minimal 
surfaces have no cusps in measure (volume density estimates); 2) minimizers satisfy the minimal 
surface equation in the viscosity sense (Euler-Lagrange equation); 3) existence of tangent cones 
(monotonicity formula). 

The main obstacles in carrying out this procedure for almost minimizers are extending the vol- 
ume density estimates and finding a substitute for the Euler-Lagrange equation. We overcome these 
first by extending De Giorgi's differential inequality argument to non-local almost minimizers, com- 
plementing the discrete iteration argument used in [5| to get uniform volume density; secondly we 
show that the almost minimality of a set forces it to solve a variational inequality, this variational 
inequality (see Section 5) is the right substitute for the minimal surface equation. 

Classically, almost minimal boundaries appear as solutions of a large class of geometric prob- 
lems: prescribed mean curvature, obstacle problems, phase transitions (for instance, their regularity 
theory plays an important role in Luckhaus' work on the Stefan problem with Gibbs-Thomson law 
[ill]). Accordingly, analogous applications arise in the non-local setting (e.g. prescribed non-local 
mean curvature). In forthcoming work, we shall use this regularity result to study the non-local 
mean curvature flow in the spirit of [l2| . The present result will also be used in the future to study 
the boundary behavior of non-local minimal surfaces. 

Remark. As s — > 1 the energy J s (see Section 2) converges to the perimeter. Thus E being 
(Ji, p, (5)-minimal means E is an almost minimal boundary in the classical sense. It is an interesting 
question whether our estimates are uniform, giving perhaps a new proof of the Almgren-De Giorgi- 
Tamanini regularity theory for almost minimal boundaries. This issue will be addressed elsewhere. 

The article is organized as follows: In sections 2 and 3 we introduce the basic definitions and 
notation, in particular the notion of almost minimal sets for non-local energy functionals, and 
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present several examples that justify this definition. In section 4 we prove the uniform volume 
estimate, which extends that of [5j to the almost minimal case. 

The key steps of the proof Theorem 1 1 . 1 1 are carried out in sections 5 and 6, where we prove the 
Eulcr-Lagrange inequalities (Theorem 15. 3|) and the Partial Harnack estimate (Theorem I5.6j) these 
imply an improvement of flatness via compactness and blow up techniques. The strategy of the 
proof is discussed at greater length at the beginning of each of these sections. 

Finally, in section 7 we prove a monotonicity formula and state without proof several propositions 
that imply by known methods the estimate on the singular set ( Theorem II .2[) . 

2. Definitions and notation 

First, we define what it means for a set to be <5-flat. 
Definition 2.1. A set E C K" is said to be 5-flat in B r (x^) if there is a unit vector e such that 

dEnB r (x ) C {\(x-x ) -e| < Sr}. 
Definition 2.2. Given any two measurable sets A,Bc M™, and s € (0, 1), we shall write 

L S (A,B):= [ [ {l-s M^^ dxdy. 

Note that the integral above might take the value +oo. It easy to see that if E C K n then 

\\ X E\\ 2 Hs/2 =c niS L s (E,E c ). 

This quantity will play the part of a perimeter functional in all that follows, hcuristically one 
can think about it as a sort of "interfacial energy" , just as the perimeter measures the amount of 
interaction between different states near the interface (see and the references therein. With 
this in mind, given a domain f2 we consider a "localized energy" functional which measures the 
contribution of fl to the H s / 2 norm of \e- 

Definition 2.3. Let n > 2 and let Q C W 1 be a bounded domain with Lipschitz boundary. Given 
E C fl, we shall denote by J s (E;fl) the quantity 

J S (E; VI) := L S (E n Q, E c ) + L S (E\ Q, E c nn). (1) 

We will from now on drop the s subscript from L s {., .) without causing too much confusion. 
We are ready to define the notion of almost minimality for the functional J s , it extends the notion 
introduced by Almgren for the perimeter (see [1|). 

Definition 2.4. Let s 6 (0,1), 6 > and p : (0,6) — > M + . We say that a measurable set E is 
( J s, p, 6) -minimal in D, if for any xq £ dE, any set F and any r such that r < min(<5, d(xo, dfl)) 
and FAE C B r (xo) we have 

J s (E;n)<J s (F;n)+p(r)r n - s . (2) 
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The assumptions we impose on the modulus of continuity p are standard, they can be found in 



Almgren's Memoir jl| p. 96a and Tamanini's monograph [15[ p. 9. They are as follows 



Al) p : (0,(5) — > R + is a non decreasing, bounded function. 
A2)p(t)=o(l). 

A3) The function t s p(t) is non increasing in (0, 5) and for some m > n + s we have 



<5 

z^ 1 (p(z))™ dz < +oo. (3) 



We emphasize that these are entirely analogous to the assumptions made for hypersurfaces of 
almost minimal perimeter. In particular, for any C > and a G (0, s] the assumptions above are 
always satisfied by p(t) = Ct a . If a p does not satisfy the first part of A3, this means pit) < Ct s for 
all small enough t, thus any ( J s , p, (5)-minimal boundary is also ( J s , Ct s . <5)-minimal, so our theory 
also applies to this casfl The following auxiliary function will play a central role in our theory. 

Remark. Throughout this work, we will implicitly assume 8 > 1 to simplify things. All argu- 
ments hold for arbitrary small S > with few modifications (and absolute constants will depend 
on 5, of course). Moreover, by rescaling one can carry the result for S > 1 to the general case. 

Definition 2.5. Given p satisfying the assumptions above, we define an auxiliary function p. 

ft 

p:(0.S)^R, p(t):=— Z - l p(z)™dz (4) 

Where m is the same constant greater than n + s which comes as part of assumption A3. 

Remarks. Note that by A3 we have p(t) — > as t — > 0, and that p is a non-decreasing continuous 
function. As we will see in Section 6, the C 1 modulus of continuity in Theorem II .11 is actually C 1 '^. 
In light of this, when p < Ct a , Theorem 11.11 savs that E has a C l a boundary in -B1/2, where 
a' = a/m. The following estimate for p will be very important at several stages of the proof. 

Lemma 2.6. Let p and p be as above. Then we have the bound 

p(te)"<m Ve£ (0,1), t G (0, S). 

Proof. By definition, 



p(te)™ mj Z p(te)-< 
Since the integrand is positive and e G (0, 1), 



p(te)~ mJo Z p(te)-< 
Thanks to condition A3, we conclude 



J a zpite)- Jo ^teJ z 



One expects that for p(t) = o(t s ) almost minimality reduces to minimality, see remark preceding Corollarv l5.4l 
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Thus, 



U^r* r 



z™~ l dz = (et)~™(te)™ = 1 



p(te)- 

Which shows that p(te)™ < p(t) for all e £ (0, 1). □ 

To finish these preliminaries, we observe that, similarly to the case of a non-local minimal 
boundary, being (J S) p, 5)-minimal is equivalent to having the following two properties: 

Supersolution property: For any xq £ dE and A C E° n B r (xo), r < min(i5, d(xo, dCt)) we have 

L(A, E) - L(A, E C \A)< p(r)r n - s . (5) 

Subsolution property: For any x$ £ 9-E and A C £H B r (xo), r < min(<5, d(xo, dfl)) we have 

L(A-B\ A) -L(A,E C ) > -p{r)r n - s . (6) 

To see the above equivalence, it is enough to check that for any set F, such that FAE CC CI, 

J a (F;tt)- J S (E;Q) = L(A-,A + ) 

+ {L(A-,E \ A-) - L{A-,E C )) - (L(A+,E) - L(A+, E c \ A+)) , 
where A + = F \ E, A~ = E \ F. 

In particular, we observe that 

J S (E; Ox) - J S (F; ^) = J S (E; fl 2 ) - J S (F; 2 ), (7) 

whenever EAF C O2 C Oi- 
Due to the this, when the set is clear from the context we may write J S (E) instead of J S (E; ft). 

3. Examples 

Boundaries of almost minimal perimeter appear as solutions of various geometrical problems 
in the calculus of variations. In this section we discuss some of the analogous problems for the 
non-local energy, our presentation extends the discussion done for the perimeter by Tamanini [l5l ] . 

Example 1. [Minimal boundaries] Trivially, any non-local minimal boundary (i.e. a mini- 
mizer of J s among all sets that are prescribed outside f2) is a ( J s , 0, 5)-minimal set. 

Example 2. [Smooth Boundaries] As it is the case for the perimeter, any set E with a 
smooth boundary is ( J s , p, <5)-minimal for p and 6 depending on the continuity of the normal to dE. 

In particular, let E C M. n be the subgraph of a smooth function u : W 1 " 1 — > M which is globally 
bounded and C then E is a ( J s , Ct a , <5)-minimal set in any compact set of K n , for some a £ (0, s). 
Of course the constants C > and 5 > depend on u and the compact set in question. 
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Let us prove this. First, note we may assume without loss of generality that the plane {x n = 0} 
is tangent to dE at 0, so that u(0) = 0, Vit(0) = 0. Let F be a measurable set such that FAE is 
contained in a ball of radius r centered at some point, which we may also assume to be the origin. 
Then, thanks to ([7]) we have that 

J s {E;n) - J s (F;n) = J s (E;B r ) - J s (F;B r ). (8) 

Regardless of what the compact fi is, as long as B r C fi. On the other hand, under the regularity 
assumption of u it is straightforward to check dE n B r C {|a; • e n \ < Cr 1+a }. In parallel to the 
perimeter, one may use a monotonicity formula (see Section 7) to show that for any F and small 
enough r we have 

Js(F;B r ) > (1 -Cr a )J a (H;B 1 )r n -", 

where H is the horizontal half-space going through the origin. Additionally, by the regularity of u 
it can be checked (with a tedious but simple computation) that 

ME; B r ) < r n - s J s (H; B x ) (1 + 0(r a )) , 

the last two inequalities imply that E is then ( J s , Ct a , 5)-minimal for some 5 > and C > 0. 

Let us sketch this last computation. It is convenient to do the change of variables x — > rx, 
E — > E r , we check easily that 

J S {E] B r ) = r n ~ s J s (E r ; B x ), 

where E r is now the subgraph of the function u r (x) — r~ 1 u(rx). To compute J s (E r \B\) we are 
going to change from the variables (x,y) to the variables (v,w), as follows: 

x = T r (v) = v + u r (v')e n , y — T r (w) = w + u r (w')e n . (9) 

In other words, T r flattens the boundary of the set E r down to the half space H = {x n < 0}. 
It can be easily checked that 

DT r = I + Vu r ® e„, det(DT r ) = 1. 

Therefore 

J s (E r ;B 1 )=L(E r r\B 1 ,E c r ) + L(E r \B 1 ,E c r nB 1 ) = 

i . . r r i 



H" JHnT, 



E? J Er n Bl \x - y\ n+s dxdy + Je^ Bi J Er \B 1 \x - y\ n+s dxdv 
1 f f 1 

- — — — — dvdw + / - — -— ~ — dvdw 

VCBi) \T r {v) - T r (w)\ n+S JH-nT^B,) JH\T r (Bi) \ T r{v) - T r (w)\ n + S 

A{v ^ r) -dvd w + I I dvdw. 



H" JHnT^Bi) \ v — w\ n+s J H c C\T r (Bi) J H\T r (Bi) \ v ~ w \ 

Where 

A(vwr) \v-w\ n+s K1 \u r (V) - u r (w')\ 

A[v,w,r)- lTr{v) _ Tr{w)ln+ s < l + L \ v _ w] ■ 



\n-\-s 
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For a C > depending on the bounds on u. From here, it is not hard to see that 



J s (E;B r ) <r n - s {l + Cr a )J s {H-B x ), 

just as we wanted. 

Example 3. [Boundaries of least energy lying above an obstacle] Let E C K™ solve 
the obstacle problem in Q with respect to a set L, i.e. it minimizes J s (.) among all sets containing 
L n ft and agreeing with a given fixed set outside ft that contains L \ O. The set L is known as the 
obstacle and the set dE n L is known as the contact set. If L is (J s ,p, <5)-minimal, then so is E. 

To see this, let F be a set such that EAF C B t (xq) C fi. Since F U L contains L and agrees 
with E outside f2 we see that F U L is admissible, therefore by the assumption on E 1 

ME) < J S (FUL). 

Now, one may check easily from the definition of J s that 

ME u L) + j s (f n l) < M F ) + M L )> 

therefore 

ME) ~ ME) < ML) - ME n L) < p{r)r n - s , 
where the last inequality follows from EAF C B r (xo) and the assumption on L. 

Remark. Later we will see that near the contact set we can always apply Theorem 11.11 an d 
thus dE is always C 1 near the obstacle L as long as L is smooth enough. 

Example 4. [Boundaries with prescribed non-local mean curvature] Another important 
class of examples of almost minimal boundaries consists of the minima of functionals of the form 
Per(.E) + J E j(x)dx. These appear for instance in phase transition problems where the mean 
curvature of the interface is related to the pressure or the temperature on the interface. It can be 
shown that the mean curvature of such minima (understood in a weak sense) is given by 7(2), when 
instead of Per(£') we consider the J S (E) energy what happens is that the non-local mean curvature 
of E must agree with 7. 

Consider the functional F(E) = J S (E) + J E ^f(x)dx with 7(2;) G LP (p > n/s). Suppose Eq is a 
minimum of J- among all sets E that agree with Eq outside J7, then Eq is ( J s , Ct a , 5)-minimal in 
fi, with a = s— ~. 

To see Eq is almost minimal, let F be another set, if EAF lies in a small enough neighborhood 
of E, then F{E) < F{F), which gives 

M E ) + / l{x)dx < M F ) + / l{x)dx 

JE J F 

J S (E)-J S (F)< / j(x)dx- / ~/(x)dx < / \j(x)\dx 
Jf J e Jeaf 

by Holder inequality, we get J S {E) - J S {F) < |£AF| 1_ p (J EAF \-f(x)\ p dx)*. 
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Therefore, if EAF c B r (xo) for some point xq and some small r, we get 

ME)- J S (F) <C nPl {r)r n ~ s 

where p 7 (r) < C n r s ~~ ||7||l*>- 

4. First estimates: almost minimal boundaries have no cusps. 

The first step in the regularity theory of minimal surfaces consists in showing that if xq G dE 
and E is minimal, then dE separates any small ball centered at xq in two regions of comparable 
volume, so that at least heuristically dE should not have cusps. In this section we prove such a 
volume estimate for non-local almost minimal boundaries. 

Proposition 4.1. Let E have the subsolution (super solution) property in f2 and xq G dE D f2, 

there are constants ro = ro(n,s,p) and cq — Ca(n,s,p) such that for any x G dE and any r < 
ram(ro,d(xo,dQ)) we have 

\EDB r {x)\ > Co\B r \, {resp. \E C n B r (x)\ > c \B r \ ). 

In particular, if E is (J s , p, S) -minimal we have 

\EDB r (x )\ 
C °- \BM\ £1 " C0 ' 

Remark. The proof is done following De Giorgi's differential inequality argument, which is 
also used to prove uniform volume estimates for minimal surfaces (cf. section 2.8 in [7[ for further 
discussion), it exploits the two inequalities available to us, namely Sobolev's embedding and the 
subsolution (supersolution) property of E. Our proof is slightly different from the one used in [B[ 
where a discrete version is used, both are nonetheless in the same spirit. 

Proof. We will prove the first of the two inequalities, the second one is obtained mutatis mutandis. 
For r G (0,1), set A r := E n B r and consider the function 

f{r) := f'(r) := H n ~ x (E n dB r ). 

The identity for /' is standard and follows from the co-area formula. By the Sobolev inequality, 

2 71 

||u||z,p<C„ >s ||u|| j, p= , 

H 1 n — s 

for u — XEnB r i we obtain 

f(r)^ <CL{A r ,A c r ). 
Since E has the subsolution property, and A r C SI n E n B r (xo), 

L(A r ,E c )<L(A r ,E\A r )+p(r)r n - s , 

so, 

L(A r , A c r ) = L(A r , E c ) + L(A r , E \ A r ) 
< 2L(A r , E\A r ) + p(r)r n - s < 2L(A r , B r (x ) c ) + p{r)r n - s 
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To estimate the L(A r , B r (xo) c ) term, we note that for all i£i r we have the bound 

r i f°° i 

/ \x v \n +s dy<C n J ^z^dz<C n Ar-\x\)- s , 

JBS. \ x y\ Jr-\x\ z 



this yields to 



L{A r ,B°) = I I X ^^ X S^ ) dxdy<C n , a I f'(z)(r - z)- s dz. 



Hence, we are led to the integro-diffcrcntial inequality, 

f(r)^ < C n , s (j* f'(z)(r - z)- s dz + p{r)r n 



Since we have /' on the right hand side, we may integrate the last inequality in the r variable 
on the interval (0, t) and get 



f(r) — dr<C n>s U I f'(z)(r~z)- s dzdr+ I p(r)r n - s dr 
interchanging the order of integration, 

ft r r 

f'(z)(r - Z y s dzdr = I f'(z) I (r - z)- s drdz, 



10 JO 

we get the integro-differential inequality, 

ft 

.1-8 



(I 



f(r) — dr<C n , s t l - s f{t)+ p(r)r n - s dr. (10) 



Let g(t) — Cot" — J * p(r)r n 1 dr, with a conveniently chosen Co(n, s). It is not hard to check 
that then we have (with the same C„ iS as in (flUl) ') 

g(r)^dr > 2C„. s t 1 - s 3 (t) + 2 f p(r)r n - s dr V t < r Q {n, s, p) (11) 

Jo 

Claim. From the integro-differential inequalities (TlOl) and (TTTI) (or equivalently ((4]) and (| 1 1 [) ) we 
conclude that f(r) > g(r) for all r such that B r (xo) C f2. 

If this claim is true then we get the lower bound on /, since due to property A3 of p there are 
c o = co(n, s, p) and rg — ro(n, s, p) such that g(r) > cor n whenever r < rg. 

Proof of the claim. We argue by contradiction as when dealing with differential inequalities: let 
r\ be the supremum of the set {r : f(t) > g(t) Vi < r}, assume that the claim was not true in 
which case r x < mm{r ,d(xo,dil)}, then for a sequence hk —> 0, hk > we have thanks to (|10p 
and CCE]) that 

g(n + h k ) > fin + h k ) > Ct s - X j /( r )(»-)/» _ p( r )r n - s dr + 2g{r x ) 

=S> liminf g(n + h k ) > 2g(n), 
which cannot be since g is continuous, and this proves the claim. □ 
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Remark 4.2. We have in fact proved something stronger than what the statement of Proposition 
I4.1l savs: under the same hypothesis as before, there is a constant tq — ro(n, s, p) such that for any 
x € dE and any r < min(ro, d(xo, dO,)) we have 

\E n B r {x)\ > c(n, s)r n - f p(z)z n - 1 dz, 

Jo 

(resp. \E C n B r (x)\ > c(n, s)r n - f p(r)z n - x dz ). 

Jo 

The volume density estimate implies (by a standard argument) that L 1 convergence of almost 
minimal boundaries {Ek} guarantees their uniform convergence, this will be useful in the last section 
when we consider blow ups and tangent cones so we state it now as a corollary. 

Corollary 4.3. Let {Ek} be a sequence of (J s , /?(.), 5) -minimal sets such that for some measurable 
set E we have 

E k ^E mLl c (R n ) 
then Ek — > E uniformly in compact subsets of il. 

We end this section with a compactness estimate , the proof is the same as that of Theorem 3.3 
in Q and it relies on the lower-semicontinuity of J s with respect to L 1 convergence. 

Proposition 4.4. Suppose that for each k > the set Ek is (J s , pk, 5) -minimal in f2, assume also 
that pk(t) — > p(t) (pointwise) and {pk} satisfies the assumptions A1-A4 uniformly in k. Then if 
{Ek} converges in L] oc (£l) to a set E, this set must be ( J s , p, 5) -minimal in Cl. 

Remark. Observe that whenever we have a sequence of almost minimal sets with pk as above 
we may use Sobolev embedding to get a subsequence that converges in L . 

5. Euler-Lagrange inequalities and partial Harnack 

In this section we will show how the almost minimality property of a set forces it to satisfy 
what we will call "Euler-Lagrange inequalities" . These play the same role that the Euler-Lagrange 
equation (in the viscosity sense) plays in the theory of minimal surfaces, as shown first in [4| for 
the perimeter and later in [5[ for non-local energies. As in the latter work, the Euler-Lagrange 
inequalities will lead us here to a Harnack estimate whenever there is enough flatness. 

We shall motivate the main idea by proving a maximum principle in a very simple case, after 
this, our main effort will be to see how close we are in the general case to this particular ideal 
situation and how can we adapt the argument below to a general setting. 

Maximum principle. Suppose that E is a minimizer of J s (.) among all sets that agree with 
E outside fi, and suppose that 

{x : x n < 0} \ n C E. 

That is, suppose that E outside VL lies on the same side of some hyperplane, then it must lie on the 
same side of that hyperplane everywhere 

{x : x n < 0} C E. 
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Proof. We argue by contradiction, to do this, let us slide from below a plane parallel to {x n = 0} 
until we touch dE inside f2. Since E does not contain {x n < 0}, there exists some c > such that 
H = {x : x n — — c] touches from below dE at some point xq € f2. Then, for all small enough e > 
consider the sets: 

Aj = E c n {x : x n < -c + e} 
A+ = E c nT € (A~) 
A e = A+UA-, 

here T e is the reflection with respect to the hyperplane {x : x n = — c + e}. The set A e will be used 
to perturb E, note that because the hyperplane is touching dE we have \A e \ > for all small e. 
Now, since A e C E c , and E is a minimizer, we may use ([5]) with p = to get 

L(^ e ,£)-L(,4 e ,£ c \A e ) <0. 

T e preserves the kernel \x — y|~("+ s ) because it is an isometry, therefore we can use it to compare 
above the opposite contributions of E and E c \A e to the quantity above and get some cancellations, 
so let F e — T e (E c \A e ). Clearly, F e C E and the inequality above is equivalent to 

L(A e , E\F e ) + [L(A e , F e ) - L(A e , T(F e ))] < 0. 

Moreover, A e can be decomposed as S e U D e , where T(S e ) = S e and D c C A~ , so we can rewrite 
things one more time: 

L(A e ,E\F e ) + [L(S e ,F e ) - L(S e , T(F e ))] + [L(D € ,F e ) - L(D e ,T(F e ))} < 0, 

by the symmetry of S e we have L(S e ,F e ) — L(S e ,T(F e )). On the other hand, because F e ,D € C 
{x n < -c + e} it can be checked that L(D €7 F € ) > L(D e ,T(F e )), thus 

L(A e , E \ F e ) < => L{A C , E \ F t ) = 0. 

Since \E \F e \ > for all small enough e, this means that \A e \ = 0, a contradiction. So it must be 
that {x : x n < 0} C E as we wanted to prove. □ 

This argument tells us that if E is a minimizer in say Bi(0) and E \ -Bi(O) is trapped between 
two hyperplanes: E C {—a < x n < —b} then the same holds inside i?i(0). It suggests we could 
do a similar argument if instead of sliding a hyperplane we could slide a very flat smooth surface 
from below (say a large ball or flat parabola); using such surfaces as barriers we may show arguing 
again as above that if E is flat enough in a cylinder, then in a smaller cylinder E cannot stretch too 
much vertically; this is part of the intuition behind the proofs of the Euler-Lagrange equation and 
the partial Harnack estimate of this section. First we recall the Euler Lagrange equation proved in 
Q (also Theorem 5.1 there): 

Theorem 5.1. Let E be a minimizer in Q, suppose that G dE and that E n contains the ball 
Bji(—Re n ), R>1. Then there exist vanishing sequences 5k, tk (with Sk << tk) and sets Ak, with 
A k C B Sk (0) n E c , such that 

L(A k ,E\B Sk ) - L(A k ,E c \B Sk ) < (C^r s + 4)\A k \, for some r, > 0. 

Moreover, 

\t \n+s — 
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As we discussed previously, the Euler-Lagrange equation implies the partial Harnack inequality 
for flat minimizers. The same perturbation argument provides a similar estimate for almost mini- 
mal sets which comes with an extra term of order p(e)e n ~ s . 

Remarks. We are going to do a slightly more general version of this estimate, which takes 
into account the deviation from minimality. Firstly, we will need to say something more about the 
sequence €k — > along which the estimates hold, in some sense, we will prove that in fact these 
inequalities hold for any e small enough (see Theorem I5.3[) . this feature will be necessary later in 
the proof of Lemma 16.41 

Another important obstacle in these kind of arguments is that the perturbation sets A t might 
have a lot of eccentricity (large diameter, but with low density). This means that the change in 
non-local energy can be several order of magnitudes larger than the available volume bounds for A e . 
In fact, A e may have a diameter of order e whereas we only know in general that \A e \ is only > e 2 ™, 
thus, one could be concerned that such perturbations do not detect the almost minimality property. 

In an earlier version of this manuscript an attempt to fix this issue consisted in constructing a 
different perturbation than in Q, however such a construction did not solve the problem either, as 
was pointed out to the authors by the anonymous referee. After this realization it was observed 
that possibly larger terms arising from the original perturbation argument can all be absorbed into 
p(.), and this "technicality" is completely taken care of by Lemma [2.61 This is still not entirely 
satisfactory, for the price we pay is not being able to prove that the first variation of J s energy 
vanishes when p = o{t s ) which is to be expected (see remark after Theorem I5.3|) . however, these 
estimates are more than enough for the regularity theory. 

Construction of the perturbation. Suppose B = B2r{— 2i?e„) is a ball contained in E 
touching dE at 0, assuming always R > 1. For any e e (0, R), we will denote by T e (x) an in- 
volution which should be thought of as a reflection with respect to a perturbation of the sphere 
dB R+£ (-Re n ). 

More precisely, T e (x) is the unique point x\ lying on the ray from —Re n to x such that the mid 
point of the segment xx\ lies on a slightly deformed sphere 3Vr j£ , where 

V RyE = {xe R" | \x + Re n \ <R + d £ (x)}, 

where d e (x) = e 2 d(x/e), d(x) = - \x'\ 2 ) + . 

Algebraically, the map is given by 

T e (x) +x = -2Re n + 2(R + d e (x)) * + ^ n . 

\X "T JXC n | 

Observe that T e is a smooth diffcomorphism and an involution of the interior of the "ring" 
V2B,e \ Vb, e into itself. We are going to perturb the set E by adding to it A e , defined by 

A e = A~ U T e {A- ) \ E where A~ =V R , e \E. 

An observation that will be important in all what follows is that A e can be decomposed as 

A e = S € UD e , 
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where T e (S e ) — S e and 

D e C V R , e \ E. 

The other properties of A e and T e that we need are summarized in the following proposition. 

Proposition 5.2. Assume that B 2 r{— 2Re n ) is an interior ball tangent to dE at with R > 1 and 
that T e is as above. Setting r x c := d(x, dVn :e ), the following holds V e £ (0, (6n) _1 ) 

1. If x e V2R, e \ V .e we have H 

|Z)T e (:r) — P x | < ~(3nr s , t + \x'\), 
where P x = reflection along the direction x + Re n . In particular, 



\T e (x)-T e (y)\ 



f - y\ 



< — max{3nr X)£ + \x'\,3nr y ^ + \y'\}, 

tc 



whenever x,y G V2_r, £ \ Vb, e . 
2. W^e /icwe i/ie inclusions A~ c i?2e B e 2 /(2R) \ E c A e c Bg e . 

Proof. To show the first assertion we are going to quantify how far is T c from being a isometry. 
First let us differentiate T e , 

x f 1 + 2 (fi+d^x)y) [I _ 2 (R + d e (x)) x + Re n ^ x + Re 



\x + Re n \ 

+e(Vd)(x/e) (g) 



\x + Re n \ \x + Re r , 
(x + Re n ) 



Re,, 



\x + Re n \ 

From this expression we observe already that when \x + Re n \ = R + d £ (x) (i.e. when x 6 9Vr, c ) 
then DT e is nothing but a reflection (specifically, P x ) plus a "defect" term Vd{x) ® if^f^r- In 
general, note that \x + Re n \ = R + d e (x) ± r x<£ and 



£>T £ (x) = P x + 2 



i? + d £ (x) ± r^e 
+e(Vd)(x/e) 



1 /-2 



a; + i?e„ x + Re r< 



\x + Re n \ \x + Re r , 



(x + Re n ) 



Since 



(R±e 



R+d e (x)±r :c 



R.+d c (x)±r :c 



\DT e (x)-P x \ < in 



\x + Re n \ 

-, and |Vd(x)| < 2R~ 1 \x'\ we obtain 



R + d e (x) -r Xjt 



2\x'\ 6n 21a;' 



R ~ R 



R 



provided e < R/2. Now, since P x is an isometry, by integrating along the segment between x and 
y we get (for all x, y in the ring under consideration) 

\T e (x)-T e (y)\ < (l + ^max{3nr x ^ + \x'l3nr yye + \y'\}^ \x - y\. 



here we are using the Euclidean norm for matrices \A\ = Tr(Aj4*) 



1/2 
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This gives the desired upper bounded for |Te( fa)l ■ Given that T e is an involution, we get 
the lower bound too and we are done proving the first assertion. 
Now to the second part: by construction, 

A' CV R>e \B 2R {-2Re n ), 

in other words, if x 6 A~ C V R . C \ B2 R (—2Re n ) then 

\x + Re n \ <R + ed(x'/e), \x + 2Re n \ > 2R. 

From these two inequalities follows that d(x' /e) ^ =>■ \x'\ < e. Moreover by a straightforward 
computation it can be checked that x n > — i? _1 |a;'| 2 , and recalling e < R we get x n > — e, the 
upper bound x n < e is obvious. All this proves that \x\ < 2e, that is A~ c i?2c- Now using the 
estimate on T e we obtain 

3n + 1 

\DT e \ < 1 + 2^— e<4. 

R 

Here we used again e < (Bn)" 1 and R > 1, this gives a bound for a ball trapping T^A^). In 
fact, we have proved that A e c i?8e- Finally, since i? > 1 we have |x| < e 2 /(2R) => d e (x') > 
3e 2 /(4i?) which implies |x + i?e n | < i? + d e (x'), this proves that Vr !£ contains B e 2/( 2 R), and thus 
-Se 2 /(2i?) \ E C an d w e are done with the proof. □ 

We are now ready to state and prove the sharpening of Theorem 15 . 1 1 that we need in our setting. 

Theorem 5.3. (Euler- Lagrange inequalities) Let E be a set satisfying the supersolution property 
f3P in with respect to (J s ,p, d), p satisfying A1-A3. Suppose that 6 dE and that En£l contains 
the ball B 2R (-2Re n ),R > 1. 

There are constants Co(n,s,p) and ro(n, s, p) such that: 

whenever < 8e* < 5 < min{ro(n, s, p),d}, 
one can find at least one e € (e*/2,e*) such that: 
L(A e , E\B S )- L(A e ,E c \ B s ) < C (p(8e)e"- s + iT 1 ^ |^|.) (12) 
An analogous inequality holds when E satisfies fifjj) and i?2i?,(^e„) C E c n fl. 

Proof. We divide the proof in 3 steps. In the first step we use directly the definition of the superso- 
lution property ([5]) to show that the left side in (fT2|) is controlled by p(.)e"~ s plus a term containing 
the singularity of the kernel (i.e. the higher order part). In the second step we show there are 
cancellations in the extra terms due to the symmetry of S e under T e that help us get rid of the high 
order singularity, this is where the built-in symmetry of A e becomes crucial. Finally, we use again 
the supersolution property to control whatever is left by L(A ( T, CBn +e (— Re n )). 

Step 1. Fix e,<5 £ (0,d) with 8e < S so that A € C Bs by Proposition [O] Define 

F £ , 5 :=T e ((E c \A e )r\B s ). 
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The first step consists in proving the inequality 

L{A e ,E \ B S ) - L(A e , E c \ B s ) < p(8e)(8e) n - s 

- [L(A e ,F ei5 ) - L(A £ , T(F ei5 ))} . (13) 

To do this, we shall use (J5J) for A = A e , with r = 8e. We decompose the integrals in ([5]) into 
parts corresponding to integration inside and outside the ball Bs, as follows 

L(A e ,E \ B S ) - L(A e ,E c \ B s ) < p(8e)(8e)"- s - I , 

Io = L(A e ,E n B S ) - L(A e , (E c \ A e ) n B s ). 

By construction, it is obvious that F €t s C E. Furthermore, we claim T e maps Bs \ Vr jC into Bs. 
To prove this let us take 

x e B s \V R , e . 

Then, we observe that the normal hyperplane to x + Re n passing through —Re n + (R+d e ) |^^"| 
separates x and the origin, so that the distance between and x is no smaller than the distance 
between and the reflection of x by the aforementioned hyperplane, this image is just T e (x). This 
shows |T e (x)| < \x\ < S and the claim is proved. We conclude that F ei s C E D Bs and 

I = L(A e , Er\B s \ F e>s ) + [L(A e , F e , B ) - L(A e , T e (F eiS ))} , 

given that the first term is non-negative, (|1 3|) is proved. 



Step 2. Using A e = S e U D e , we rewrite the integrals on the right hand side of (fTl?|) , 

L(A e ,F €tS )-L(A e ,T e (F e , 5 )) = [L(S e , F e , tf ) - L(5 6 ,T e (F e , tf ))] 

+ [L(D e ,F etS )-L{D e ,T e (F e , s ))]. 

Let us bound each term. We apply the change of variables x, y — > T e (x),T e (y) to estimate the 
integral on S e x T e (F Et s) as follows 

= f I ( lct(DT.(x)),et(DT(y)) 

Jf c , s Js< \T e {x)-T e (y)\«+' 

The first part of Proposition [5T^1 asserts that \DT t {x) — P x \ < 2i? _1 (inr x ^ + \x'\) for all x under 
consideration, which gives an upper bouncjfl for the Jacobians, leading to 

L(s e ,T e (F ei5 )) <(i-s) j Ft 5 j Se 1+c -*~ lm g;;j^ 

Also by Proposition 15.21 we may bound \T e (x) — T e (y)\~ n ~ s from above with \x — y\~ n ~ s , times 
a factor given by Proposition [521 namely 



\T e (x) - T e (y)\- n ~ s < [1 - 2R- 1 max^m^ + 2\x'\,3nr y , e + 2|y'|}]-"- s |x - y\ 



3 under our assumptions \DT t {x) -I\<2, thus \dct(DT e )\ < (1 + \DT e {x) - P x \) n < 1 + n2 n \DT e (x) - P x \ 
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Putting these two bounds together we conclude that 

L(S £ ,T e (F £ , s ))<L(S £ ,F e ,s) 

+(1 - s) L A — ^ dxdy - 

For the second term we also do a change of variables but only in y, using again Proposition 15. 2 
to bound the Jacobian. However, we control the denominator in a different manner, arguing as 
in the proof of the claim in step 1 to conclude \x — T e {y)\ > \x — y\ for all x, y 6 Vr i£ (remember 
x € D e ,y £ F e ,s Vr_ c ). In this case, 

L(D e ,T e (F e<5 )) = (l-s) [ [ 1 dxdy < 
<(!-*) / / 

JF, * J L 



1 + C n i? _1 max{3nr x e + 2\x'\, 3nr v e + 2\y'\} , , 

— — dxdy . 



iF e , s Jd c \x - y\ n+s 

The bounds we have obtained for L(S e ,T e (F t ^)) and L(D t ,T e (F e ^ j) give us, 
L(A e ,F et5 )-L(A e ,T e {F et s)) > 

-C 1 - s ) If.,, U ° nR " m ^ 3 -;r+^> 3 " r "' +2|y ' l} dxdy. 

Now, denote by the set A t \ A~ , observe that T e (A+) C A~. It is reasonable to try to 
control L(Af, F e j) in terms of L(A~ , F Cl s), that is, using the change of variables x — > T e (x), 

L { At^ s) ^l- s) J F J AtJ -^dxdy 

f f l + C n R~ 1 m^{3nr x . e + 2\x'\,3nr y . £ + 2\y'\} 

< (1 - s) / / : — |— r- dxdy, 

Jf, . J at I* - 2/|" +;s 



and we get 



L{A e ,F etS )-L(A e ,T e (F eiS )) > 

-(1 - fl) J F / A _ g»«" 1 °'"{3»r, 1 «+2|x'|,3nr H ,«+2| 1/ '|} da;dy _ 



Note that the quantity in the numerator is of order at most S (since the integration is over 
subsets of Bs), which is not sharp enough for our purposes. Nevertheless, when x and y are far 
apart the kernel is bounded so in that region we can get a better bound. Specifically, we have 

to , oi 'i o , oi H^ / f 6n F - y| + 4e if \x — y\ > 2r x , e 
max{3mv e + 2M,3nr,, £ + 2M}< j ^ . f _ Jj < ^ . 

This follows at once from the inequalities 

r y , e <r x , e + \x-y\, \y'\<W\ + \x-y\ 
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=> 3nr y . e + 2\y'\ < 3nr Xie + 2\x'\ + An\x — y\, 

Together with the fact that r X}<L < e and \x\ < 2e for all x £ A~, which was shown in Proposition 
[5721 Then it follows that 

r max{3nr x , g +2|:r'|,3nrv g +2|t>'|} , 
JF t ,s \x~y\ n + s & 

< C n (fF e:S \B 2rxe (x) \x-y\™ + ° d y + /f.^HB^ jt (x) \x-y\™ + ° d v) • 

Integrating in x G A~ (and noting that \x — y\- r <- s + 1 is locally integrable), this gives 

, f f max{3nr x , e + 2|a; / |,3n.ry i e + 2|y'|} 
(1 - s) / / !-j ,— — cft/dir. 

<C n (5 1 - s \A:\ + eL{A-,F e , 5 )), 
and using our previous bound for Jo we get 

-Jo < CnR- 1 (lA^lS 1 -* + eL(A-,F e:S )) . (14) 

This finishes the second step, next we show that L(A~ , F £i s) is not too big. 

Step 3. Recall from Step 1 that F e j C E. Therefore L(A~ , F e< s) < L(A~,E), and this last 
quantity (by the supersolution assumption is not bigger than 

L(A;,E c \A;)+p(8e)(8e) n - s . 

Moreover, since E c \ A~ C e we conclude 

L(Aj,F 6tS ) < L(AJ, V£ e ) + p(8e)(8e)"- s . 

CZazm 0: There exists eo(n, s, p) £ (0, 1) such that for any e* < eo there exists some e s (e*, 2e*) 
such that 

Again following [5j, we (essentially) rewrite L(A~ , e ) in terms of 

o(r) := H"" 1 (A7 n 9Vr ;7 .) , < r < e. 
Indeed, if a; € A~ C Vr^ then for d{x) := d(x, dV^ e ) 

1 f°° 1 

V c^ \x - y d(x) r"+ s 

1 , „ »I-Bi| 



/ \ X -y\n+ S dy - S( l ir y 



4 Compare with Lemma 5.2 in 
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Thanks to this, we get the bound 

L(A 7 , V£ e ) < C n>a [ 1 dx, 

since d(A~ , <9Vr j£ ) < e and d(x, 9Vr, £ ) ~ (e — r) if x E 9VR, r (r < e) this implies that 



Jo \ e r ) 



(e - r) s 

Then, let us suppose that the Claim is not true, then there exists a small e* such that the 
following inequality holds for every e E (e*, 2e*) (for brevity we write rj = ^-) 

(e-r)' 



-<ir > e / a(r)dr. 



We integrate this for e E (e*, 2e*) and change the order of integration to get for f(r) — | Vr.,- \ E\ 
defined for r E (0, e) (recall that f'(r) = a(r) by the coarea formula) 



2e* 

/'(r)(2e* - rf- s dr + I f'{r) ({2a - r) 1 - 3 -{a- r) 1 " 5 ) dr 



1 - s 

> 



f'(r) [(26*) 1 -" -r 1 -"] dr+ 



i-n 

+ I f(r) ((2a) 1 -" - a 1 "'') dr 



Observe that the second term on the left hand side is bounded from above by 2 1-s /(e*)(e*) 1 ~ s 
and the second term on the right hand side of the inequality is bounded from below by C s f(e*)(e*) 1 ~ n , 
this shows that if e* is smaller than some eo = £o(s) then 

26 f(r){2e*-r) X - s dr > C s J ' /'(r) [^e*) 1 -" - r 1 ""] dr 
=► (e*) 1 - 3 /'(r)dr > C^CO 1- " / ' /'(r)dr. 



Observe that we can use the same proof of Proposition 14.11 to get density bounds for E c with 
respect to the modified balls Vr, € . In particular f(r) satisfies an integro-differential inequality like 
(|10p. from where it follows easily that J„ e f'(r)dr / for all positive and small enough e*. This 
means we can cancel the common factors and get (e*)^ -5 > C s . Recalling -q — -Mp > s we get a 
contradiction if e* is too small and the claim is proved. 

Using the bound for L(A~ , F £t $) together with step 2 we obtain the bound 

-h < CnR- 1 {s 1 - s \A-\+e^ ± \A-\+ep(8e)e n - 

for some e E (e*, 2e*) as long as e* < eo, and as long as S < ro for some ro = ro(n, s, p). Moreover, 
since e < 5 we will for the sake of brevity bound the term e^ - ^^ by (J - 2— | | . Putting all this 
together with the bound (fT3|) from step 1 the theorem is finally proved. □ 
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As a corollary of this proof, we get a pointwise Euler-Lagrange inequality whenever p(t) < t p 
(p > n + s). This is far from a sharp estimate, as it should hold for any p > s, however, this 
requires a more refined perturbation argument than that of Theorem 15.31 The main justification 
in presenting this corollary is that it will introduce an estimate used later in the proof p. 

Corollary 5.4. Suppose E satisfies the supersolution property |3P in Q with respect to a p such 
that p(t) < Ct p ,p > n + s. Then anytime E has an interior tangent ball at a point xq 6 Q we have 

XE (x) - xe : (x) dx < 
\x-x \ n+s 

As in the previous theorem, we get the opposite inequality in the case of a subsolution and an 
exterior tangent ball. 

Remark. Given xq g dE, the quantity 

-{is) I x ^ x) - XE : {x) dx. 

'J Rn \x-x a \ n+s 

is what is known as the non-local mean curvature of dE at xq. It has many features similar 
to the usual mean curvature, for instance, thanks to the invariance of the kernel \x — y\~ n ~ s under 
rotations, it is easy to see that a ball of radius r has constant non-local mean curvature and that 
this constant is proportional to r~ s . Further, bounded convex sets E would have boundaries with 
positive non-local mean curvature. 

Proof. Without loss of generality, we may take xq — 0. Let e and 5 satisfy the hypothesis of 
Theorem 15. 31 observe that 

L(A ei E\B 5 )-L(A e ,E c \B s ) f X e(x) - X e°(x) 



C n e 



< / rr^—dx < C n e5~ 



-l-s 



This follows from the Lipschitz bound for \x — y\ n s in each of the sets {|a; — y\ > 8} and the 
fact that A t C B% t . Then Theorem 1 5 . 31 implies that 



Recall that by Propositions 14.11 and 15.21 we have \A~\ > Coe 2 ™, therefore 

I XE{X Lir iX) ^ < C n eS~^ + C i (Come— + 6*) 



5 a posteriori one might use the knowledge that dE is smooth to prove this for all p > s using for instance 
perturbations by smooth vector fields. 
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since /(8e) < C n e p , p > n + s we are free to let e — > with 6 > fixed, getting 




dx < Cq—S 



l-s 



taking the limit as 5 — > the corollary is proved. 



□ 



To finish the section, we prove the Harnack estimate. This is the key tool needed for the 
improvement of flatness argument. We will use the following definition (see also discussion at the 
beginning of Section 6): 

Definition 5.5. A set E C R™ is said to be flat of order a at x with respect to 7 : [0, f] — > R + , if 

a = j(2~ k ) for some k £ N and there is a sequence of unit vectors {e{\i such that 



Remark. As we are assuming that p is defined in (0,6) with S > 1 (see remark before Definition 
2.5) we may take the auxiliary function p (again see Definition 2.5) as our 7 above. 

For the purposes of the next section (namely showing regularity) it will also be convenient to 
state the Harnack estimate in a different scale, i.e. scaling the ball of radius B 2 -k(0) into the 
ball i?i(0), if a set E is ( J s , p, d)-minimal, then the rescaled set will be ( J s -pk, rffc)-rmnimal, for 
p k {t) = p{2- k t) and d k = 2 k d. 

Theorem 5.6. (Partial Harnack) For any p satisfying assumptions A1-A3, d > and s G (0, 1) 

there are positive constants ko(n, s, p, d) and So{n, s, p, d) with the following property: 

Suppose that E is (J s , p k , dk)-minimal in B 2 k (0), flat of order a = p(2~ k ) with respect to p and 
that G dE. Then, if k > fcg the following two inclusions hold: 



Here Pk(t) — k t),dk '■= 2 k d and eo comes from the flatness hypothesis. 

Remark. This theorem says that conditioned to dE being "fiat enough" near zero one has a 
Harnack inequality. Note that unlike the usual Harnack inequality for regular elliptic equations, we 
cannot reapply it over and over at finer scales (since the flatness condition might fail in a smaller 
ball), thus the term "Partial Harnack". We get back to this point in Section 6. 

Proof of Theorem ] 5. 6{ Suppose that for given k, a — p(2~ k ) and 8q there is a set E satisfying the 
assumptions and such that one of the two inclusions does not hold. We may assume that eo is the 
positive direction on the x n -axis. With this setup, we shall get a contradiction by showing that if 
a and S are picked universally small then E will "stretch" too much, forcing it to contradict the 
Euler-Lagrange inequalities (Theorem 15. 3|) . 

Without loss of generality (the other case is dealt with similarly), suppose that 



dEr\B 2 i(x ) C {\(x- x ) -ei\ <2 l ~/(2 l - k )}, VI: 0<l<k. 



E n B So C{x-e < o(l - <$o)}> B to n {x ■ e < -o(l - 6%)} C E. 



B s n {x n < -a(l - 5 2 



(15) 
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Slide from below a ball B of radius R=a~ 1 > 1 until it becomes tangent to dE from the interior, 
which we can do thanks to the flatness assumption. Since we are assuming (fT5|) . without loss of 
generality the ball must be tangent at least at one point z such that 

z e B s n {x n = -a(l - S 2 )}. 

We are under the assumptions of Theorem 15 .31 thus using the same S as in (|15[) and any e so 
that e and S are allowed, we have 

L(A e ,E \ B 5 {z)) - L(A e ,E c \ B 5 {z)) < C (p(2- fc 8e)e"- s + aS^ \A~ |) , 

our goal is to obtain an opposite bound. Note that B C E so for any y G K n we have 



Xe(x) - Xec\aS x ) 



Bx\B s {z) 



\x-y 



n+s 



dx > 



Xb{x) - Xbc\aS x ) 



i\B s (z) 

2" 



\x-y 



n+s 



dx 



Xe\b(x) 

\n+s ' 



dx. 



Bx\B s (z) \X-V 



Since G dE and z G B$ n {x„ < — a(l — 5 2 )}, the ball B a+2 s(y) contains E f) B a any time 
|z — y| < (5. This means that the function \x — y|~("+ s ) must be greater than c n (S + a)^^ n+s ^ in 
E n B&, which is a subset of (E\ B) D (Bi \ Bg(z)) if a and S are small enough and a ~ <5 (recall 
z lies on {x„ = — a(l — S 2 )} so Bg(z) and £? are away from the origin). Thanks to Proposition ^. II 
this set has measure no less than cpa™ when a is small. Hence we have 



rr _ „ n+s — "V / 2 w 



Additionally, we can check by direct computation that 

Xb(x) - XBc(x) 



Bi\B s ( z ) 



1 3 - y| n+s 



da; > -CJ- S Vy G 



Then, since Xe(x) ~ Xe c \a € ( x ) — Xb — Xs c we conclude there is a universal Co such that 
whenever C^ 1 a < S < Co a we have the inequality 



Xe(x) - xe-{x) 



bi\b s (z) \x-y 



n+s 



dx>Co L S~ s \JyeB L {z). 



(16) 



Next, we estimate this same integral outside Bx (0). Let \y\ < 1/4 and r > with 2'~ 1 < r < 2 l 
for some I such that < I < k. Since E is flat of order a with respect to p, it is not hard to prove 
that 



Xe{x) - Xec(x) 



x—y\—r 



\x-y 



n-\-s 



dS 



< 



aC(n, p) 



for some 770 = r]o(n, s, p), this is in part due to assumption A3. We integrate this with respect to r, 



Xe(x) - Xe-(x) 



2 



\x-y 



n+s 



dx 



< 



aC(n, p) 



»t+rio 



dr + C„ 



„l+s 



dr 
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< (J Cor-^dr) a+^2- sk , V y e Bi. 



Since k > fco, choosing fco large enough (equivalently, choosing a small enoug h) to get 2- sk < 
a = p(2~ k °) we obtain 



Xe{x) - Xe4 x ) 



Bi (0) c 
3 



F - y 



< C a, Vy G Bi 



(17) 



The estimates (fTrJj) and (|T7)) put together give us (as long as 8 « 1) 

f Xe(x) - xe- (x) 

Jb 5 (zY \x-y 



-dx > 8 s — Co(n, s, p)a, V y s Bs(z). 



From the way we picked e we can make sure that A c C Bs , this means that we can integrate 
the above inequality for y € A e and reach the lower bound: 

L(A € ,E \ B s (z e )) - L(A tl E c \ B 5 (z e )) > 8- s \A e \ - C{n, s, p)a\A e \. 



It is worth emphasizing that this last estimate is independent of the Euler Lagrange inequalities, 
it follows only from the flatness hypothesis, the density estimate from Proposition 14.11 and the 
assumption that (Ti~5)) does not hold. Putting these two bounds together we see that for all e, 8 and 
a under consideration we have 

8~ s \A e \ -C(n,s,p)a\A e \ < C Q (p(2- fc 8e)e"- s + S 1+1 ^ \A^\\ . 



We now apply Lemma [2701 which says that if 8e < 1 then 

p(2- k 8e) < p(2- k ) m = a m , 

so if a, e and 8 are universally small (e being admissible for Theorem I5.3[) . 

<T S < C [a m e n -'|4 s |- 1 + l] . 

To finish, recall that |A e | > cot 2 ™ (Propositions 14. ll and [5^]) and note we can finally pick all 
parameters so that e ~ 8 ~ a, we get 



since m > n + s, this means that 



<T S < C* [S m - n - s + l] , 



8- s < C . 



Choosing 8 (universally) small we get a contradiction. This implies the estimate since we have 
shown that the first inclusion (|15l) must hold with some (universal) 8 when fco is (universally) large, 
since fco grows as a — ¥ 0. The other inclusion is dealt with in the same way, except we must slide a 
ball B from above and use the subsolution version of the Euler-Lagrange inequality. □ 
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6. Flat boundaries are C 1 



The concept of improvement of flatness is nowadays well understood, it arises not only in the 
regularity theory of minimal surfaces but on free boundary problems as well. To achieve it here, we 
make use of a method developed by Savin to address the regularity of level sets in phase transitions 
See also [l3j where a similar idea is applied to elliptic equations. 

Heuristically it goes as follows: by a standard argument, if fiat boundaries are not C 1 then 
there is a sequence of (J s ,p,d) minimal sets Ek, with vanishing flatness in B 1 / 2 k, and such that 
dEk cannot be trapped inside a flatter cylinder in B x / 2 k+i. For this sequence there is a partial Har- 
nack estimate (in our case Theorem l5.6[) which, as we will see below, implies that if we dilate each 
surface in the direction in which they become flat, then a subsequence of them converges uniformly 
to the graph of a continuous function u with controlled growth at infinity. This function is then 
shown to solve the equation (— A)^~u — in all of K™ which shows it must be linear. This must 
give a contradiction, since the uniform convergence and smoothness of the limit force the sequence 
to lie eventually inside flatter and flatter cylinders, against our initial assumption. Therefore the 
original sequence cannot exist. 

There are new technical problems in carrying out this procedure for non-local minimizers, and 
they were dealt with in |5| . The main difficulty comes from the influence of far away terms in all 
the estimates. This would not be a problem if we knew that the surface was flat enough away 
from a neighborhood of 0, but after scaling a much smaller neighborhood of into one of size of 
order i we arc sending new parts of the surface far away and the non-local terms could become 
relevant again. In 5] it is shown this is not really a problem, by taking into account not only a 
small neighborhood of but also larger and larger balls where the flatness grows until it becomes 
of order one. This is actually the same as the standard flatness hypothesis, except we are keeping 
track of the information between the scale of order 1 and the smaller scale. Later in the argument, 
as we rescale, this information allows us to control the non-local terms. This idea was already used 
in the statement and proof of the Partial Harnack estimate (see also Definition 5.6). 

In concrete terms, the main result of this section, which implies Theorem 11.11 as in is the 
following: 

Theorem 6.1. Let E be (J r s , p,d) -minimal in B±, £ dE Pi B\ and p satisfying A1-A3. There is 
a universal ko = ko(n, s, p) such that if for some k > ko we have a sequence of inclusions 

dEf]B 2 -i c {\x-ei | < 2- l p[2- 1 )} VI < k 

for some sequence of unit vectors {ei} l<k , then we can find unit vectors e; for all I > k such that 
the inclusions above remain valid for all I. Recall p is the auxiliary function defined in Section 2. 

To understand the asymptotic behavior of these sets in a small ball, we are going to change 
the scale by a factor of 2 k to make the ball B 2 -k into the unit ball. Then the previous theorem is 
equivalent to: 

Theorem 6.1. (Rescaled) Let E be (J s , pk, dk)-minimal in B 2 k with pk(r) = p(2~ k r), dk = 2 k d, 
G dE with p satisfying A1-A3. There is a universal ko = ko(n, s, p) such that if for some k > kg, 
the set E is a-flat at with respect to p, a = p(2~ k ), then we can find a unit vector e_i such that 
we have the inclusion 
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8EnB 1/2 C {\x-e_ 1 \<2- 1 p{2- 1 - k )} 



(18) 



The proof of the theorem will be divided in a couple of lemmas, all of which deal with a sequence 
of sets {Ek}. Here for each k > the set Ek is ( J s , p/e, d^) -minimal as above and each being dfe-flat 
with respect to p, — p{2~ k ). 

(k) 

In order to normalize things further we also assume for all k that agrees with the unit vector 
in the positive direction of the x ra -axis (we can always get in this situation via a rotation). 

To the sequence {Ek}k we can apply directly the partial Harnack estimate from Theorem 15. 6[ 
as the classical oscillation lemma, it gives a Holder estimate: 

Corollary 6.2. There is a universal kg such that if k > ko, then for any x G dEk, the set 
dEk H B$(x) can be trapped in between the graphs 

{(y', Vn) ■ Vn = x n ± Ca k max{a^, \y' - x'| 7 }}, (19) 
where C > and 7 £ (0,1) are universal constants. 

Proof. The proof is exactly the same as the proof of Holder continuity for solutions of elliptic 
equations via the Harnack inequality. That is, if v solves an elliptic equation in B\ then applying 
Harnack inequality k times we see that 

sup v — inf v < p I sup v — inf v 1 , p, € (0, 1). 

B 1/2 k B l/2k \B 1/2 k-l B l/2k-l J 

Since k is arbitrary, one concludes the oscillation of the function decays geometrically and thus 
v is C 7 at with 7 = log 2 p. The difference is that in our case we can only apply the Harnack 
estimate only as long as the flatness assumption holds (thus it is a "partial" Harnack estimate). 
More specifically, we can apply partial Harnack for a fc-th time as long as ap k < ao, where p = 1 — 5$ 
(Sq as in Theorem l5.6|) . thus the maximum number of times we can do it for the set Ek is ~ log(|^) 
and from here the Holder estimate follows. 

□ 

This corollary is the key in proving the following lemma, in which we show that from any 
sequence of vertical dilations of Ek , one may always pick a subsequence converging to the graph of 
continuous function, for a function that does not grow too much at infinity. 

Lemma 6.3. For each k we dilate the set Ek vertically, and define 

El = {(x',x n ) : (x',a^ 1 x n ) G E k }, 
then, considering the new sequence {E^}k, we have: 

1. Along a subsequence {E^k converges uniformly in compact sets ofM™ to the subgraph of a 
Holder continuous function u : R™ -1 — > R 

2. Moreover, it(0) =0 and for a universal C we have 

\u{x')\ < C(l + \x'\p(\x'\)), V x' G R"- 1 . 
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Proof. Since the functions \E k are bounded in H s ^ 2 there is a subsequence that converges locally 
in L 1 to some set. By the partial Holder estimate, and the fact that a k — > 0, this convergence is 
actually uniform in compact sets. Since E dE k for all k the same is true for the limit set. Again 
by the partial Holder estimate, we see that the limit set can be touched above and below at by 
the graphs {x n — ±C|x'| 7 }, its not hard to see after a translation that we may do the same at every 
point of the boundary of the limit set. Therefore this set is the subgraph of a Holder continuous 
function u : R' i_1 -> E. 

Now we prove the second statement, u(0) = is clear since E dE. Moreover, by a diagonal- 
ization argument we may assume that for each fixed I that the sequence of vectors {e\ } converges 
to a unit vector e/. If we denote by pi E R™ -1 the horizontal projection of e/, we can see that 

\u{x') -pi-x'\< C2 l p(2 l ), if x' E B 2 , , 

and \pi-pi+i\ < C2 l p{2 1 ), V I > 0. 
Therefore, \u(x')\ < C(l + \x'\p(\x'\)), W £ which finishes the proof. □ 

The last lemma we need says basically that when we "linearize" the non-local minimal surface 
operator, we obtain the fractional Laplacian, and thus, the limits of non-homogeneous blow ups 
must be harmonic. 

Lemma 6.4. Let u be the same function from the previous lemma, then 

(-A)^u = 0«il n - 1 . (20) 

Remark: Given the growth condition of u at infinity, it follows by estimates for elliptic equations 
that u is linear (see for instance Landkof's treatise [10(). 

Proof. We shall prove that (|2Tfl) holds in the viscosity sense (cf. Q). Let ^ be a smooth function, 
and to fix ideas, suppose it is touching u from below at the origin. Fix positive numbers ei, M and 
K; Lemma [63] says that one can find some k > K such that dEk fl Bm lies in a a^ei-neighborhood 
of the graph of auu{x), ak = p(2~ h ). Moreover, it is easy to see that dEk must be touched from 
below by a vertical translation of the graph of ak<p{x) at some point Xk E B ei (0), in particular 
dEk is being touched at that same point by a ball of radius Rk > C^a^ 1 . 

3 /2 3/2 

For each k let us pick so that a k < £fc < 2a k and such that €k is admissible for Theorem 
15.31 We will approximate 

L(A eh , E k \ B Sk Or)) - L(A eh , E c k \ B Sk (x)) 
with the fractional Laplacian of u. In fact, we make the following claim. 

Claim: 

JB> H (v)\B' t {V) \x '-x '\ n + s -2a k J BS \x-x\"+° 

+ 0{af) + 0(M-" 2 ) + C(0)<5 1 - s , for universal Tj% : T)% > 0. (21) 
Recalling the proof of Corollary 15 Al we have that 

XE k (x)-XE t (x) < £j _ (1+s) L{A t ,E k \B s )-L{A t ,Ej\B s ) 
fl . (s) \x-x\ n +* ~ " \A t \ 
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Theorem 15.31 guarantees that for all large enough k and fixed small 5, 

XE k ( x ) -XE?(x) 



B k (x k ) \X-X k \ n +* 



-dx < a k C Q (a^ 1 e fc J- (1+s) + C^ 1 ' 



+ (a^p(2- k 8e)e- s ) =: a k C B k . 

In light of this and (|2~T|) to prove the lemma it suffices to get a bound on B k ■ Observe that from 
the way we picked e k the first term of B k obviously goes to zero as k —> oo as long as 5 remains 
fixed, since 

-1 ^ o -1 3 / 2 i) 1/2 v n 

a k e k < 2a k a k ' = 2a k ' -> 0, 

the second term in B k remains bounded in k so we only need to focus on the third one. Once again 
we apply Lemma 12.61 which says that if k is large then 

al 1 p{2- k 8e)el s = p(2- k 8e)a k 3s/2 < p(2- k ) m p(2- k )- 3s / 2 , 

by property A3 (cf. Section 2 ) we always have m > 2 and we conclude 

lim a k - l p(2- k 8e)e k s = 0. 

k— >oo 

In other words, for each fixed 6 > we have lim Bk < CqC^S 1-3 . Putting this together with 

k— >oo 

TJ) we get for 8 > fixed, 

U !- X ?~^ dx' < CoC^-* + o(l)) + 0(ei) + O(M^). 

I B' M (x' k )\B' s (x' k ) \ x X k \ 

Taking t\ — > 0, fc — s> +00, we see x' k —> 0, and the inequality above gives us 



L 



u{x,) ~ U ^dx'<C5^ + 0(M-^). 



B' u {0)\B' s {0) W\ n+S 

Then letting M — » +00 and 5 — > we obtain the desired inequality. The only thing left to prove is 
the claim (|21l) (which was already proved in Q). 

Proof of the claim. Consider for any 5 > the cylinder 

D m (x) = {\x' — x'\ < m; \x n - (x) n \ < m} C K™. 
Using the fact that in D c & (x) the function \x — x\~ n ~ s is Lipschitz with constant C(6), 

1 1 

\x - x\ n+s \x' - x'\ ri 

Recalling also that dE k lies in a a^ei-neighborhood of the graph of a k u(x), we get 



< C(5)a 2 , whenever \x n — (x) n \ < 2a. 



[ XE j X) -* E J {x) dx = 2 [ a fc Kx0-,(xQ + O( £l ))^ + 

Id m (x)\D s (x) \x~x\ n + s J B ' M (x)\B' 5 (x) \X'-X'\ n + S 
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This is almost estimate (|2"Tj) . except we still need to take into account the integration over 
and D$ D Bg, . We use the fact that dEk is being touched from below by the graph of ak<fi, thus we 
get the lower bound 



XE k (x) - XE? (x) 



\n+s 



dx > -a fe C(0)5 1 



I D { (x)\B e (x) \X~X\ 

Next, note we can neglect the contributions outside Dm{x) 



XE k {x) -XE?(x) 



\n+s 



dx 



< 



a k r n -"p{r)r- n - s dr + Ca 



fc 



< Ca k M- r ' 2 + Ca k 



and with this we are done. 



□ 



Proof of Theorem \6.1\ If the theorem were false, there would be a sequence of sets {Ek}k>o satis- 
fying the hypothesis of the theorem but such that for each k, the set dEk n -B1/2 does not lie inside 
any cylinder of flatness 2~ 1 /3(2~ fe ~ 1 ). To this sequence we readily apply Lemma [6.31 and Lemma 
16. A\ so passing to a subsequence we see that the dilated sets E% converge to a hyperplane. Since 
the convergence is uniform, for any e > we can find a fci for which dE% n B\ lies inside an e 
neighborhood of the limiting plane. 

Taking e small enough and going back to the set E^, we see that this set contradicts the 
assumption that none of the sets dEkV\B X j 2 could be trapped inside cylinder of flatness 2~ 1 /3(2 -fe ~ 1 ), 
and this proves the theorem. □ 



7. Monotonicity formula and regular points 

In the classical theory of minimal surfaces, an important role is played by the monotonicity 
formula, which says that whenever E is a minimal surface: 

lf^)>0, ifOEE. 
dr \ r"- 1 / ~ 

In this section we consider an analogous formula for the non-local case, introduced in We 
discuss briefly how it is used to characterize the regular points of the boundary of an almost minimal 
set E, several proofs are omitted. 

Definition 7.1. For any function u £ L 1 (M",^) where /1 = (l + |.t| 2 ) 2 dx, we define the 
extension of u, u : -> E, as the solution to 

f div(z a Vu) =0 in R'| +1 

\ u = u ond = W l 

Here R™ +1 denotes the half space {X = (x,z) : x £ R", z > 0}, and a = 1 — s. 
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For more about this extension and its connection to integro-cliffcrcntial equations, see [2]. In [5] 
it is shown how the non-local energy J S (E;Q) relates to the z a dxdx) energy of u when 

U = XE - XE-- 

Proposition 7.2. (See section 7.2 in (H) Assume E and F are two sets such that J S {E; f2), J S {F; ft) < 
oo, and suppose that EAF CC O. Then 

inf f z a (\Vv\ 2 - \Vu\ 2 ) dxdz = c n , s (J S (F;Q) - J S (E;Q)) . 
F Jn+ 

We will omit the proof of this proposition, since there is nothing to be added to the proof in 0| 
for it to work in our situation. With this fact in hand we can prove the monotonicity formula for 
almost minimal sets. 

Lemma 7.3. (Monotonicity formula) Suppose E is (J s , p, d) -minimal in tt, and that G dE. Let 
Ue = Xe — Xe c , then the function 

L+z a \Vu E \ 2 dX r 
®E(r) = JBr ' + (n-s) p(t)t n - s - l dt (22) 



is monotone increasing for all small enough r. Moreover, this function is constant if and only if E 
is a cone with vertex at and p = 0. 

Proof. We can compute $' E (r) in the standard way to get 

&e(t) = 4~s{ I z a \Vu-E\ 2 dX - / z a \Vu E \ 2 dx) + (n - s)p(r)r n — J . 

r \JdB+ r JB+ J 

The theorem will follow at once from 

I z a \Vu E \ 2 dX> [ z a \Vu E \ 2 dX - p(r)r n - s . (23) 

n - s JdB+ JB+ 

To prove (f2"3"]l . consider the function 

v(X) - / ^(rX) if |X| < r 
{ ' \ u(X) if \X\>r ' 

where we are putting X — X/\X\. Note that the trace off) on {z = 0} equals xf(x) —Xf^(x) where 
F is some set which agrees with E outside the ball B' r (0) C M™. Then, by Proposition 17.21 and the 
almost minimality of E we have 

/ z a \Vv\ 2 dX> [ z a \Vu\ 2 dX - p{r)r n - s . (24) 

JB+ JB+ 

Using the construction of v, we compute: 

z a \Vv\ 2 dX= f f (tz) a \(Vv)(tX)\ 2 t n dS dt 

B+ JO JdB+ 
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VupQ = r\X\^ [l-X®Xj (Vu E )(rX) = r\X\- 1 {V s u E )(rX), 
in the last identity Vs is the component of the gradient tangent to the sphere. Then, 

/ z a \Vv\ 2 dX = r 2 [ t n+a - 2 [ z Q |(V s u E )(rX)\ 2 dS dt 

JB+ JO JdB + 

t n+a - 2 dt z a \{V s u E ){rX)\ 2 dS = - / z a \{V s u E )(rX)\ 2 dS. 



/as+ n — s ./3R+ 

Changing variables in the last integral, we get 



L 



z a \\7v\ 2 dX = — ^— f z a \(V s u E )(X)\ 2 dS. 
n-s J dB + 



This together with (|24p implies (j2"3"]l and the monotonicity is proved. Actually, we have proved 
the sharper estimate 



*i(r) > -L (f Wdx) 



From here it also follows that if $>' E (r) = then u u = on every dB r , which shows that u 
must be homogeneous of degree 0, in particular it must be so when restricted to {z — 0}, but there 
u — xe ~ Xe c , so E must be a cone with vertex at the origin. In this case by the scaling of the 
energy we see that p = 0. □ 

With the monotonicity formula in hand, one can reproduce easily the blow up analysis of classical 
minimal surfaces (see Q for the classical case). The main steps are summarized in the following 
theorem whose proof we omit. It is worth mentioning that here one needs the results at the end of 
section 4. 

Theorem 7.4. Let E be a (J s ,p,d) -minimal in B\, and G dE. 

1. Suppose rk — > is a sequence such that 

E rk -)• C in L\ oc 

then C is a minimal cone with vertex at 0. Moreover, the blow up sequence E rk converges 
locally uniformly to the minimal cone C. 

2. We have an energy gap: there exists a constant Sq such that if C is a minimal cone which 
is not a half-space, then <I>c > &h + $o, where H is any half-space. 

3. If the minimal cone C is a half-space, then dE is C ' p in some neighborhood ofO. In particular, 
if Sq is as in (2) and $^(0) < <&h + <^o then dE is smooth in a neighborhood of 0. 

Finally, let us mention that the dimension reduction analysis done in [5( directly applies to our 
situation. Thus the estimate achieved there for the dimension of the singular case carries over to 
almost minimal sets too. Therefore, thanks to Theorem ll.il we have: 

Theorem. If E is ( J s , p,d) -minimal in VL, then the Hausdorff dimension of the singular set T, E C 
dE fl Q is at most n — 2. 
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Additionally, the monotonicity formula and the theorem above show that, in the case of the 
obstacle problem (see Section 3), all points of E near the contact set are regular, specifically, we 
have the following result: 

Theorem. Let L CC Q be a domain with a smooth boundary. Suppose E minimizes J s (.; 0) among 
all sets containing L, then dE is C 1,a in a neighborhood of L for some a < s. 
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